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Using explicit solutions of the QCD transport equations, we derive an effective potential for an 
anisotropic quark-gluon plasma which under plausible assumptions holds beyond the Hard Loop 
I—} , approximation. The configurations, which are unstable in the linear response approach, are Char- 

on ■ acterized by a negative quadratic term of the effective potential. The signs of higher order terms 

I ' can be either negative or positive, depending on the parton momentum distribution. In the case of 

a Gaussian momentum distribution, the potential is negative and unbound from below. Therefore, 
(f) ■ the modes, which are unstable for gauge fields of small amplitude, remain unstable for arbitrary 

^ | large amplitudes. We also present an example of a momentum distribution which gives a negative 

■ quadratic term of the effective potential but the whole potential has a minimum and it grows for 

' sufficiently large gauge fields. Then, the system is weakly unstable. The character of the instability 

is important for the dynamical evolution of the plasma system. 

o . 
in ■ 
o ; 

^ : 

Ph| I. INTRODUCTION 

An anisotropic quark-gluon plasma is unstable with respect to the transverse mode known as the filamentation or 
^3 Weibel instability [1-6]. The instabilities have recently focused some attention in the context of the phenomenology 
| of nucleus- nucleus collisions at high energies. The experimental data on the particle spectra and the so-called elliptic 
• *~J ■ flow, which have been obtained at the Relativistic Heavy-Ion Collider (RHIC) in Brookhaven National Laboratory, 
suggest, when analyzed within the hydrodynamical model, that an equilibration time of the parton 1 system produced 
at the collision early stage is as short as 0.6 fm/c [7]. Calculations, which assume that the parton-parton collisions 
are responsible for the equilibration of the weakly interacting plasma, provide a significantly longer time of at least 
2.6 fm/c [8]. To thermalize the system one needs either a few hard collisions of the momentum transfer of order of the 
characteristic parton momentum 2 , which we denote here as T (as the temperature of equilibrium system), or many 
collisions of smaller transfer. As discussed in e.g. [9], the time scale of the collisional equilibration is of order 

* hard ~ g 4 Hl/g)T ' (1) 

where g is the QCD coupling constant. The characteristic time of instability growth is roughly of order 1/gT for a 
sufficiently anisotropic momentum distribution [1,4,6,10]. Therefore, the instabilities are much 'faster' than the hard 
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x The term 'parton' is used to denote a fermionic (quark) or bosonic (gluon) excitation of the quark-gluon plasma. 
2 Although we consider anisotropic systems, the characteristic momentum in all directions is assumed to be of the same order. 
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collisions in the weak coupling regime. Since the system's momentum distribution becomes more isotropic when the 
unstable modes develop, the instabilities have been argued [1,10] to effectively speed up the thermalization processes. 
Very recent numerical simulation [11], where the fields, particles and color charges are treated classically, fully confirms 
the argument. 

We note that the isotropization should be distinguished from the equilibration process [1]. The instabilities driven 
isotropization is a mean-field reversible phenomenon which is not associated with entropy production. Therefore, the 
collisions, which are responsible for the dissipation, are needed to reach the equilibrium state of maximal entropy. 
The instabilities contribute to the equilibration indirectly, reducing relative parton momenta and increasing the 
collision rate. And recently, it has been observed that the hydrodynamic collective behavior, which is evident in the 
experimental data [7], does not actually require local thermodynamic equilibrium but a merely isotropic momentum 
distribution of liquid components [10]. 

The stability analysis of the plasma system is usually performed within a linear response approach which assumes 
smallness of the gauge field amplitudes. An elegant qualitative argument [12] suggests that non-Abelian non-linearities 
do not stabilize the unstable modes as the system spontaneously chooses an Abelian configuration in the course of 
the instability development. This happens because the Abelian configuration corresponds to the steepest decrease 
of the effective Hard Loop potential. However, the phenomenon of abelianization which emerges from the numerical 
simulations [11-15] is more complex. In 1 + 1 dimensions, i.e. when the gauge potentials depend on time and one 
space variable, the simulations [12,13] of the Hard Loop anisotropic dynamics [16] and the classical simulation [11] 
show that initially the abelianization works well but it is less efficient when the fields become non-perturbatively large 
in the course of instability growth. The latter effect is nicely seen in the calculations [13]. In the simulations [14,15] of 
the Hard Loop dynamics [16] performed in 1 + 3 dimensions, the abelianization stops working at the non-perturbative 
scale and the growth of field amplitudes changes from exponential to linear. 

The question arises what happens beyond the Hard Loop approximation, whether the higher-order non-linearites 
stop the instability growth, or maybe they make the system even more unstable. In attempt to address the question, 
we derive an effective potential of the plasma configuration which is known to be unstable at the Hard Loop level. We 
follow the method developed in our earlier paper [17] which uses explicit solutions of the QCD transport equations. 
To get the solutions, the system under study is assumed to be static and translation invariant in two space directions. 
The effective potential, which is obtained, thus corresponds to the simulations in 1 + 1 dimensions. It holds at any 
order of the gauge coupling but a kind of Abelian approximation, which is justified due to the effect of abelianization, 
is needed to derive it. The higher order terms of the effective potential, which are relevant at sufficiently long times 
of instability development, carry information about the actual character of a plasma configuration. Depending of the 
momentum distribution of partons, the plasma is either weakly or strongly unstable. 

Although we do not even attempt to solve the dynamical problem of temporal evolution of the unstable plasma 
- we merely derive the effective potential - we estimate various time scales of the instability growth. In particular, 
we show that the quartic and higher order terms of the potential become important at the same time scale when the 
fields starts to significantly influence the particle's momenta, as explained in the concluding section. 

Throughout the article we use the natural units with Ti = c = ks = 1- In Section II and in Appendix A, where 
we follow the four-dimensional notation, we use the metric convention (1, —1, —1, —1) and we distinguish lower and 
upper Lorentz indices. In the remaining parts of the paper, where we mostly deal with three-vectors, all vectors are 
contravariant and to simplify the notation we do not longer distinguish lower and upper indices. 



The distribution functions of quarks (Q), antiquarks (Q), and gluons (G) are assumed to satisfy the collisionless 
transport equations 



II. THE METHOD 




(2b) 



(2a) 



(2c) 



where {..., ...} denotes the anticommutator; the covariant derivatives and act as 



0^=8^- ig[A^{x), ... ] , 



V tl =d f i- ig[Af,(x), ... ] , 



2 



and A^ being four-potentials in the fundamental and adjoint representations, respectively, 

A»(x) = A^x)r a , A^(x) = -if abc A»{x) , 

and fabc are the structure constants of the SU(JV C ) group. Since the generators of SU(A r c ) in the adjoint represen- 
tation are given by (T a ) bc = —ifabc, one can also write A^ = A£T a . The field strength tensor in the fundamental 
representation is — d \A V — d v A^ — ig\A^ A„], while T^y denotes the strength tensor in the adjoint representation. 

As already mentioned, the instability of interest is a very fast phenomenon. The characteristic time of instability 
development is not only much shorter than the characteristic time of hard parton-parton collisions (1) but it is also 
shorter than that of soft collisions which control the dissipation of the color degrees of freedom. In equilibrium, the 
collisions with momentum transfer of order g 2 T occur at the time scale [9] 

tsott ~ g*Hi/9)T • (3) 

Since the characteristic time of instability growth, which is roughly of order 1/gT for sufficiently anisotropic momentum 
distribution [1,4,6,10], is much shorter than t so f t , the collision terms of transport equations (2) can be neglected. 

If a solution of the transport equations is known, one finds the associated color current which in the fundamental 
representation is 

fix) = -| J dP P »[Q(p,x) - Q(p,x) - ±-Tr[Q(p,x) - Q(p,x)] + 2iT a f abc G bc (p,x)] , (4) 
where the momentum measure 

takes into account the mass-shell condition po = E p = |p|. The adjoint currents equal j%(x) = 2Tr[r a j M (a;)]. Since 
the current is given by the relation 

where S e s = / d 4 x£ c s is the effective action to be added to the Yang-Mills one, one obtains S e s integrating Eq. (5). 

Finding exact solutions of the transport equations (2) is in general a difficult task. However, it is possible to get 
such solutions under some restrictive conditions. Here we consider a system where both the gauge field and the 
distribution functions are invariant under some space-time translation(s), i.e. 

d ai A»(x)=0, n = 0,1,2,3, (6) 

and 

d ai Q{p,x) = d az Q(p,x) = d ai G(p,x) = , (7) 

for a fixed a,, which can involve more than one Lorentz index. If on = the system is static while for on = 1,2,3 the 
system is homogeneous - the gauge field and the distribution functions depend only on time. 
We look for the solutions of the transport equations (2) in the form: 

Q(p, x) = f(p ai - gA ai = £ { -^~A ai (x) A a2 (x) ■ ■ ■ A an (x) - ° J[PaJ , (8a) 

n _Q n - Op ai Op a2 . . . Op an 



oo 



Q(P, x) = f( Pai + gA a% (x)) = J2 (x) A a2 (x)--- A an (x) tt^tt , (8b) 



n=0 



n! 1 2 " dp ai dp a2 . . . dp ar 

isY A ^ A ^ A ^ d n f g (p ai ) 



G(p, x) = f g (p ai - gA at (*)) = £ { -^-A ai (x) A a2 (x)-- A an (x) q^q^q^ , (8c) 

where it is understood that a sum is taken over the repeated indices. The functions /, / and f g are, in principle, 
arbitrary but they can be fixed by additional considerations. The solutions of the form (8) do not assume smallness 
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of gA provided the series (8) are not terminated at a finite number of terms. We stress that all components of the 
four-vector p M are treated here as independent variables. 

One finds that the expressions (8) exactly solve Eqs. (2) when 

[D^A ai ,A aj ] =0 , = 0,1,2,3, (9) 

and 

[V li A ai ,A aj ]=0 ^ = 0,1,2,3. (10) 

We note that Eq. (10) is automatically satisfied if Eq. (9) holds. The condition (9) is satisfied when, in particular, A£ 
has only non- vanishing components in the Cartan subalgebra of S\J(N C ) which for SU(3) are in the a — 3, 8 directions 
of the color space. If A£ has only one color component, the condition (9) is also trivially satisfied. 

Having solutions of the transport equations, one finds the associated color current. Inserting the expressions (8) 
into Eq. (4), one gets 

# (*) = E Mi— A «\ (*) • • • A ^ (*) / dp ( n ) 

x |Tr[r a r Cl ■■■T Cn ] 

The current (11) provides, after integrating Eq. (5), the effective lagrangian of the form [17] 

C ^ = E ^TT)! { TrKl {x) ■ ■ ■ Aa ^ {x)] J dPp ' 

+ Tv[A ai (x)--- A an+1 (x)} f dP p^ ^M^il X . (12) 

The distribution functions /, / and f g include two helicity states of, respectively, quarks, antiquarks and gluons. 
Various quark flavour states are also included in / and /. We note that Eq. (12) with the equilibrium distribution 
functions provides in the static case (o^ = 0) the full one-loop effective potential [17]. 

To get the complete action, the lagrangian (12) should be supplemented by the gauge field 'kinetic' term — \F^ v F a ^ 
which involves, due to its non-Abclian nature, not only quadratic but also cubic and quartic terms of the potential 
A^. Since the lagrangian (12) has been derived under the assumption that the gauge fields obey the condition (9), 
the same condition must be applied to the 'kinetic' term, which then becomes effectively Abelian. 

As seen, the effective lagrangian (12) is local. In general, this as a consequence of the space-time homogeneity (6,7) 
combined with the commutation condition (9). However, it should be stressed that, as explicitly shown in [17], the 
g term of the lagrangian (12) does not require the condition (9) which is needed for the higher order terms. The 
quadratic term and its relationship to the potential obtained within the linear response method is further discussed 
in Appendix A. 

The Hard Loop action can be found for both the equilibrium [18-21] and anisotropic [16] plasmas, modifying the 
g term of the lagrangian (12) to comply with gauge invariance. Presumably, the higher order terms of the effective 
potential (12) can be also generalized in the same way, once the space-time invariance is abandoned. However, we 
are not going to follow this direction. Our purpose here is to understand the possible role of the higher order terms 
of Eq. (12), whose coefficients are uniquely fixed by powers of g and by the parton distribution functions. 

III. EFFECTIVE POTENTIAL FOR UNSTABLE CONFIGURATIONS 

In this section we discuss configurations which are known to be unstable in the linear response approach. Specifically, 
we consider the filamcntation instability which occurs in the homogeneous plasma when the momentum distribution 
is anisotropic. When the instability develops, the kinetic energy of particles is converted into field energy. The wave 
vector of the unstable mode is perpendicular to the direction with a 'surplus' of the momentum and it points towards 
the direction with the momentum 'deficit'. We choose the wave vector along the z axis and we expect to observe an 
instability when the momentum distribution is elongated in the x or y direction when compared to the momentum in 
the z direction. In Sec. II we have carefully distinguished the upper and lower indices of four-vectors. From now on, 



d n f( Pa% ) 
dp ai ■■■ dp a 



(-1) 



n+l 



d n f( Pa% ) 
dp ai ■■■ dp 



+ Tr[T„T Cl • • • T C J 



dp ai ■ ■ ■ dp a 



d n f( Pai ) + , 1)n+1 0"/(PaJ 

dp a2 ■ ■ ■ dp an+1 dp ai ■ ■ ■ dp a 
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all four-vectors are contravariant and the upper and lower indices are no longer distinguished (except the Yang-Mills 
lagrangian written in a relativistically covariant form). 

Let us now consider a geometry of the filamentation instability referring to classical electrodynamics [2] . Due to the 
'surplus' of momentum along the x or y axis, the fluctuations of current flowing in the x—y plane are amplified. Thus, 
we expect non- vanishing x or/and y current components. When the wave vector of a standing (time independent) 
wave points in the z direction (k = (0,0, k)) and the current (j) flows in the x — y plane, there are non- vanishing x 
and y components of the vector potential (A) and the magnetic field (B). Indeed, one checks that 

j = (a x k 2 sm(kz), a y k 2 sm(kz), 0) , A= (a x sin(kz), a y sin(kz), 0) , B = (— a y k cos(kz), a x k cos(A:z), 0) , 

which depend only on z, satisfy the electromagnetic equations 

VxB=j, VxA = B. 

As seen, the dynamics of the system, which is static and invariant with respect to the translations in x and y 
directions, can be still nontrivial. Therefore, we consider such a quark-gluon system. The solution of the transport 
equation is then of the form f(p — gA (z),p x — gA x (z),p y — gA y (zj). Since the filamentation mode is transverse and 
the magnetic field is responsible for the instability development, we can choose the gauge condition A — 0. Finally, 
the distribution function is assumed to obey the mirror symmetry 

/(P0,P)=/(P0,-P). (13) 
Then, the formula (12) for on = x,y gives the effective potential V e s = —£ e s as 

™ + & ^ + <(7^ + <S35P£» + ° (f/) ' (14) 

where 

(-■) = JdPp°--f{p). (15) 
If the quantity, which is averaged over four- momenta, is independent of po, the integration over p is performed and 

<->-/^-/<p). 

The function /(p) should be understood as f(E p , p). We note that only the quark contribution is taken into account 
in Eq. (14). Including the contributions of antiquarks and gluons is straightforward. We also observe that due to the 
mirror symmetry condition (13), there are only terms of even powers of A in Eq. (14). 

To study the system's stability, the effective potential generated by the quasiparticles should be supplemented by 
a potential- like contribution coming from the Yang-Mills 'kinetic' contribution to the lagrangian —\F^ v F a ^ v . Since 
we consider static configurations in the A = gauge, the lagrangian equals 

\pa rpa ni> — B a B a 

4 ^ u 2 
where the chromomagnetic field B a is defined as 

B a = V x A a + lf abc A b x A c . 
V 

The Yang-Mills lagrangian is built of the gauge potentials and their gradients. We identify as a contribution to the 
effective potential the term of the Yang-Mills lagrangian where gradients of A a are absent. This is the only term 
which survives in a homogeneous system, and it reads 

Vym = |.9 2 r fc 7 a<ie (A h A' i )(A c A e ) . (16) 
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Because of this term, the system has been argued [12] to spontaneously choose an Abelian configuration in the course 
of the instability development. We will assume that this is the case, thus justifying the Abelian-like conditions (6) 
under which the effective Lagrangian (12) was obtained. As noted in the Introduction, the simulations in 1 + 1 
dimensions [11-13], which are relevant for our study, confirm effectiveness of the abelianization. We return to this 
problem in the concluding section. 

The negative sign of the quadratic term of the potential (14) usually signals the instability. Let us consider the 
issue. Within the linear response analysis, unstable modes are found, see e.g. [22], as solutions of the dispersion 
equation 

dct[k 2 (P J - k l W - u 2 - IF(w, k)] = . 

The polarization tensor of anisotropic plasma, which was derived in [22], is given in Appendix A by Eq. (A4). 

Further, we consider a simplified situation where the wave vector points in the z direction (k = (0,0, kj) and the 
chromoelectromagnetic field has only one non- vanishing component along the x axis. Then, as discussed in e.g. [3], 
the dispersion equation for the mirror symmetric momentum distribution simplifies to 

k 2 -lu 2 -n xx (u,k) = o. (17) 

The unstable modes of interest correspond to the solutions of Eq. (17) with positive imaginary frequency. The 
existence of unstable modes is controlled by the so-called Penrose criterion [23] , which was applied to the simplified 
equation (17) in [1] while a general configuration was studied in [6]. According to the Penrose criterion, there are 
unstable solutions of Eq. (17) if there are such wave vectors k that 

k 2 -U xx (oj = 0,fc) < . (18) 

Since U xx (uj = 0, k), as given by Eq. (A4), is independent of k, the Penrose criterion (18) reads 



/ 



dPp x 



df _ Px df 
dp x Pz dp z 



> . (19) 



The respective quadratic term of the effective potential (14) is negative when 



/ 



df 

dP Px -^->0. (20) 

OPx 



The criteria (19, 20) appear to be fully equivalent to each other. Since the system under consideration is static and 
translationally invariant along the x and y directions, the solutions of the transport equations cannot explicitly depend 
on p z . Therefore, the last term in Eq. (19) vanishes identically. 

Essential information contained in the potential (14) is carried by the sign of each term. As discussed above, the 
negative sign of the quadratic term signals instability of small amplitude modes. Depending on signs of higher order 
terms, one deals with strongly unstable or weakly unstable configurations. In the next sections we consider examples 
of momentum distributions which illustrate the two situations. 

IV. STRONGLY UNSTABLE CONFIGURATION 

We consider here a specific Gaussian form of the distribution function /. Namely, the function is chosen as 



f{Px, Py ,Po) = 2 3 7r 3 / 2 y / /3(/3 - a x )(0 - a y ) pcxp(a xP 2 x + a y p 2 y - /?p§) , (21) 
where p is the quark density (p = (1)) while the parameters a x , a y and (0 > , > a x , > a y ) control the 



momentum distribution. One observes that taking into account the mass-shell constraint p — E p = ^Jp 2 + Py +p 2 , 
the distribution function equals 

f(p x , Pv ,E p ) = 2 3 7r 3 / 2 v //?(/3 - a x ){0 - a y ) pexp( - (0 - a x )p 2 x - (0 - a y )p 2 y - 0p 2 z ) , (22) 
and one immediately finds 

(„ 2 ) = t ( v 2 ) = t ( v 2 ) = JL 

2{0-a x ) 2(0 -a y ) 20 



G 



With the distribution function (21), the momentum derivatives in Eq. (14) are trivially computed and the effective 
potential reads 




The terms of higher powers of g, which are suppressed in the small coupling regime, can be easily calculated as 
well. The moments of the Gaussian distribution function can be computed analytically and expressed through 
hypergcomctric functions but the results are of rather complicated form, and thus they are not very useful. We 
observe that the moments of the Gaussian distribution function, which are present in Eq. (23), are all positive. 
Therefore, the potential (23) is negative when a x and a y are positive. We note that the higher order terms are also 
negative for a x > and a y > 0. Then, the effective potential is negative and unbound from below for an arbitrary 
magnitude of the gauge potential. We also observe that the infinite series with the first few terms given by Eq. (23) 
is convergent for any value of the gauge potential as the Gaussian function is regular and its derivatives of any order 
are finite. 

Obviously, the gauge field amplitude cannot roll down to infinity, as suggests the potential (23), and there are 
several dynamical processes that can stop the field growth, rebuilding the effective potential. We keep in mind here 
the non-Abelian non-linearities, isotropization, equilibration which operate at different time scales briefly discussed in 
the concluding section. Quantitative analysis of these effects requires a dynamic treatment which includes the back 
reaction and self-consistent generation of fields. We only mention here that if no other processes stop the instability 
development, the inter-parton collisions will certainly do the job. The collisions, which are neglected in our analysis, 
slowly drive the system towards equilibrium. And when the thermodynamic equilibrium is reached, the system is no 
longer unstable - the effective potential becomes positive with a minimum corresponding the equilibrium configuration. 

Below we consider two special cases of the distribution function (21) which are of interest in the context of phe- 
nomenology of heavy-ion collisions. When two relativistic nuclei of equal mass collide, the parton momentum distri- 
bution in the center of mass frame is initially elongated along the beam axis - it corresponds to a prolate ellipsoid. 
With the time passing, the parton momentum distribution in the local rest frame evolves due to the parton free 
streaming towards an oblate shape with the momentum squeezed along the beam direction. The two configurations 
are discussed below. 



A. Prolate shape 



When a.y = and an = a x , the momentum distribution (22) is 



f(p x , Py , E p ) = 2 3 7T 3 / 2 /3 v //3 - a\\ pcxp( - (J3 - a„)^ - 0{p 2 y +p 2 z )) , (24) 
where, as previously, a\\ < and > 0, and 

For qt|| > the distribution (24) corresponds to the prolate ellipsoid elongated along the x axis. The effective potential 
(23) for the distribution (24) equals 

One observes that when the momentum ellipsoid is elongated in the x direction (a\\ > 0), the potential (25) is negative 
and unbound from below. 
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B. Oblate shape 



When a± = a x = a y , the momentum distribution (22) is 

f(p x ,p y ,E p ) = 2V/ 2 V/3(/3 - a±) pexp(- (/3 - a ± )p\- (3pf) , (26) 
where a± < (3, /3 > and p\ = p x + p y . One observes that 

<* ■ <* " 5«> = • «> - 9 ' 

For a± > the distribution (26) corresponds to the oblate ellipsoid squeezed along the z axis. The effective potential 
(23) for the momentum distribution (26) equals 



v« = -i s V(^)'&[4 + 4] < 27 > 

When the momentum ellipsoid is squeezed in the z direction (a± > 0), the potential (27) is negative and unbound 
from below - all terms are negative. 

V. WEAKLY UNSTABLE CONFIGURATION 

The effective potential crucially depends on the parton momentum distribution. The Gaussian distribution discussed 
in the previous section provides an effective potential where the quartic and higher order terms are negative once 
the quadratic term is negative. In such a case, a mode, which is unstable in the linear response approach, remains 
unstable for gauge fields of larger amplitude. However, it is not difficult to invent a momentum distribution which 
gives a negative quadratic term of the effective potential and a positive quartic term. Then, an unstable mode in the 
linear response approximation is stabilized when the field's amplitude becomes sufficiently large. 

As previously, we consider here a static system which is homogeneous along the x and y directions and we choose 
the momentum distribution of the form 

2-7T 2 /? 2 

f(p*,Py,Po) = Y^Jp ^ V p[S( Px -V) + S( Px +V)] e"^° , (28) 

with p > and V > 0. As in the case of the Gaussian distribution, the normalizing constant is chosen in such a way 
that the quark density p equals (1). The computation of the constant is explained in Appendix B. The distribution 
(28) describes two counter-streaming beams of quarks with the momenta p x = V and p x = —V. The momentum in 

the y — z plane is distributed according to <~ exp^ — P^Jv 2 + p y + p£j . 

The effective potential (14) computed with the distribution (28) equals 

V eB = -p0 1T ^Tr[A x ] + -p0 Y+p~V ' 

----- ^ p/3 2 "- 1 2n ~ 1 Z^ P Tr[A 2 "] + • • • (29) 



(2n)\ r ' 1 + pV 

The calculation details are given in Appendix B. As seen, for 1/(3 < V < 3//3 the quadratic term is negative but 
the quartic and higher term are positive. Then, the mode, which is unstable in the linear response approximation, is 
stabilized when the field's amplitude becomes so large that the quartic term matters. 
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FIG. 1. The effective potential (31) for 5 values of the parameter c. The most upper line corresponds to c = 0.5, the lower 
one to c = 2.0, etc. 

Interestingly, the series (29) can be summed up and the result is 

VcS = W+PP) ^ t ~ PV ^ {f3gAx) ~ ^ + P9A x sh(PgA x ) - ch(PgA x ) + l] . (30) 

To further analyze the effective potential (30), we define the dimensionless potential U = (3(1 + (3V)V e g/p , the 
dimensionless parameter c = (3V > and the dimensionless variable s = (3gA x . The gauge potential A x is treated as 
a scalar, not matrix, quantity. Then, 

U(s) = (1 +c)(l -chs) + sshs . (31) 

One easily finds 

if^s 2 + ^s 4 for ±>1, 
^ e s for s > 1 . 

As seen, the potential U(s) always grows for sufficiently large s. Since the derivative of U(s) equals U'(s) = — cshs + 
schs, the potential U(s) monotonously grows for c < 1 as the equation s — cths has no solution. For c > 1 the 
potential has a minimum. For c only slightly exceeding unity i.e. such that < c— 1 <C 1, the potential is minimized 
for s m i n ~ •v/3(c— 1). For c> 1 the minimum is at s min c. The potential (31) is shown in Fig. 1. 





VI. DISCUSSION 



The effective potentials analyzed in the previous sections have been derived, assuming that the commutation 
condition (9) is satisfied. The assumption is discussed here in more depth. We also consider the hierarchy of time 
scales of the dynamical evolution of the unstable system. The discussion reveals relevance of our results. Let us start 
with the commutation condition (9). 

It has been argued in [12] that for the unstable configuration the color direction of the steepest descent of the 
effective potential corresponds to the Abelian configuration where the Yang-Mills term (16) vanishes. Therefore, 
the system becomes Abelian in the course of instability development. The suggestion has been confirmed by the 
numerical simulations [11-13] in 1 + 1 dimensions which are relevant for our effective potential. (As mentioned in the 
Introduction, the abelianization in 1 + 3 dimensions [14,15] is operative only for sufficiently small gauge fields.) Thus, 
there is a good reason to expect that the commutation condition (9) is satisfied. However, in the analysis [12], only the 
Yang-Mills (16) and quadratic term have been taken into account. So, one wonders whether the argument still holds 
when the higher order terms are taken into account. We note that if the strongly unstable configuration becomes 
Abelian at the Hard Loop (g 2 ) level, it remains Abelian for arbitrary large gauge field amplitudes. There seems to 
be no reason to change the runaway direction in the color space which corresponds to the Abelian configuration. We 
obviously assume here that initially the gauge fields are small. Thus, we believe that the condition (9) is justified for 
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the strongly unstable configurations, as that one corresponding to the Gaussian momentum distribution. When we 
deal with the weakly unstable configuration, as that one discussed in Sec. V, the commutation condition (9) is no 
longer justified when the effective potential starts to climb up. Then, the system's color configuration can become 
non-Abelian again. 

Our main objective was to derive the higher order terms of the effective potential. Let us give a crude estimate of 
the time scale when the quartic term becomes relevant. For the purpose of order of magnitude estimates, we treat 
the gauge potential as a scalar quantity A and we write down the effective potential as 

V eS ~ -p 2 A 2 + \A A , (32) 

where p 2 > 0. One immediately observes that when A 2 <~ p 2 /|A|, the quartic contribution to the effective potential 
is equal to the quadratic contribution. One estimates p and A as 

2 2 P x 4 P 

P ~9 7f , a- 5 ^3, 

where p is the parton's density and T is, as previously, the parton's characteristic momentum or energy. Thus, the 
quartic term is important when A ~ T / g. At this scale the system's energy, which is stored in the fields, is comparable 
to the total energy of particle excitations. When A ~ T / g not only the quartic term but the terms of any order are 
important and the system's dynamics becomes nonperturbative. 

And how much time is needed for the amplitude of an unstable mode to become so large that the quartic contribution 
to the effective potential is comparable to the quadratic one? When the instability develops, the gauge potential grows 
as A(t) = ae 7 *, where a is the initial value of A and 7 is the instability growth rate. Substituting the exponential 
dependence of A into Eq. (32), one finds that the quartic contribution equals the quadratic one after the time 



^-higher ^ hi 



1, J t 



7 "V|A|a 2 , 

To get a 2 , which we treat as a fluctuating gauge potential in a neutral background, we first estimate the fluctuating 
current j which was studied in [2]. One finds that j 2 <~ g 2 p 2 /N, where N is a number of partons in the volume which 
matters for the unstable mode generation. It is identified with the cubic wavelength of the unstable mode. Thus, 
N <~ p/k 3 , where k is the characteristic wave vector of the unstable mode, and j 2 ~ g 2 pk 3 . Physically, the current 
appears due to the charge fluctuations in the volume of interest. Further, the gauge potential squared is estimated as 
a 2 <~ j 2 /k 4 <~ g 2 p/k. Since k 2 ~ p 2 [1,3,4,6], one gets a 2 ~ gy/pT [10]. Because the instability growth rate 7 is also 
of order of p [1,3,4,6], we finally get 




(o 6 J ■ 



If p ~ T 3 , as in the equilibrium, the estimate simplifies to 



thigher ~ -jp hi(l/fir) • 

Let us now consider the characteristic time of back reaction (iback) which is defined as a time interval after which 
the change of parton's momentum (Ap) due to the Lorentz force (F) is comparable to the parton's momentum itself 
(T). One computes 

/•tback /"tback „jL„ 

Ap ~ / F(t)dt ~gk A(t)dt - ^— (e 7 * b » k - l) . 

Jo Jo 1 

Neglecting the unity in the last expression and using the above estimates of k, a and 7, the equation Ap — T provides 

i Ft t t 3 \ 

iback ~ ~\ — In I —s- <~ ^higher • (33) 

.9 V p V -<?V 

As seen, the time scale of back reaction and that of relevance of higher order terms of the effective potential are the 
same. This is not an accident - the two effects are actually the same. We could define the scale of back reaction time 
as the time when the parton distribution function which, according to Eqs. (8), is a function of p± g A, significantly 
depends on gA. And then, it is evident that <back = ^higher- 
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The shortest time scale of interest is that of instability growth and consequently, the early stage of instability 
development can be studied with the frozen momentum distribution. At the time scale of back reaction the influence 
of the mean field on the parton's motion cannot be longer neglected. At this scale the anisotropic system is expected 
to evolve fast to isotropy. The time scale of inter-parton collisions controls the local equilibration of the system. 
However, the scales of both hard and soft collisions, given by, respectively, Eq. (1) and Eq. (3), are much longer than 
iback (33). A plasma system produced in relativistic nucleus-nucleus collisions is subject to expansion due to the 
system's finite size, an effect that we have not considered. It should be noted that the expansion provides another 
relevant time scale in the problem. As shown in [3], it is so fast that the momentum distribution evolution caused by 
the expansion can stop the instability growth. 

We conclude our considerations as follows. Configurations of an anisotropic quark-gluon plasma, which are unstable 
within the linear response approach, can be strongly or weakly unstable, depending on the sign of higher order terms of 
the effective potential. The sign in turn depends on the parton momentum distribution. Examples of both situations 
have been given. As argued above, the non-Abclian non-linearities do not stop the instability development if the 
momentum configuration is strongly unstable. To quantitatively understand the system's dynamics at longer times, 
one has to go beyond the Hard Loop approximation, which effectively amounts to consider the effect of back reaction 
of the mean fields on the parton's motion. At even longer times, the effect of inter-parton collisions, which drive the 
system towards local thermal equilibrium, has to be taken into account. Studies of finite systems require a proper 
treatment of the expansion which may affect seriously the instability development. 
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APPENDIX A 

Let us discuss here the quadratic term of the lagrangian (12). The term equals 

C^{x) = ~K'»'A-{x)At{x), (Al) 



2 

where 



= — f dPp» . (A2) 

Only the quark contribution is taken into account here. The term (Al) should be compared to the lagrangian derived 
in the linear response analysis which is 

£lr(*) =-\j d 4 y A* (x)Il^(x - y)A a v {y) , (A3) 

where the Fourier transformed polarization tensor II'"' (k) of anisotropic plasma equals [22] 

W"(k) = 9 - [dPp^ 
W 2 J dp\ 



p a k a + i0+ 



(A4) 



As in the case of W , only the quark contribution has been written down. We note that IP" (A;) = W^(k) and 
kf J/ Tl^ v (k) = 0. We also observe that the Hard Loop action was found for both the equilibrium [18-21] and anisotropic 
[16] plasmas, modifying the action (A3) to comply with a gauge invariance. 
One finds that the coefficient (A2) can be written down as 

= W(u),ls. = 0) , 

where k = (w,k). It is worth noting that 

lim W v (oj, k = 0) ^ lim W(uj = 0, k) . 
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We also observe that IP 1 ' (a;, k = 0) is independent of u>. 

There is a subtle point concerning the calculations of n Ml/ and of /C". Because of its transversality, the polarization 
tensor (A4) can be computed in two different but equivalent ways: either the distribution function f(p) is treated as a 
function of all four components of the four- momentum p or we substitute into Eq. (A4) the function /(p) = f(E p , p) 
which depends only on the three-momentum. When K,^ is computed, the two methods are not equivalent to each 
other. However, the derivation of the effective potential (12) shows that f(p) has to be treated here as a function of 
all four independent components of p. 



APPENDIX B 



We present here some details of the calculations which lead to the effective potential (29) of weakly unstable system. 
The normalization constant C of the distribution (28) is determined by the integral 

P = J 0f26(po)6(p 2 )po.f(p x , Pv ,Po) = C J ^LJL[5(p x -V) + 6(p x +V)] , 

where E p = ^Jp x + p y + p\. It can be easily computed, using the variables (p x ,E,4>) with (j) being the azimuthal 
angle. Then, one finds 

P = J d Px [5( Px -T) + S( Px + V)} j | dE p E p e~ pE ' = e~ fir {\ + pp) . 

To get the effective potential (29), we first compute 

/ Px ® \ f d^p Px $ f d^p Px ® 

\padp x / J (2tt) 3 p dp x J (2ir) 3 E p dp x 

It should be stressed that the derivative with respect to p x acts only on the explicit dependence of f(p x ,p y ,E p ) on 
p x . Therefore, after the partial integration the derivative acts on p x /E p and on the E p dependence of f(p x ,p y ,E p ). 
Then, one finds 

V, a v pfiel" r , „ m d ( .„.-, A- 13V 



^Pod Px / 1 + PVJ ' dp x V* > " 1 + pP 

It is not difficult to generalize the above result to the higher order terms. Keeping in mind that 

dPx 

one finds for I = 1, 3, 5 . . . the following result 



^Pod Px / 1 + pV 

which gives the effective potential (29). 
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